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Abstract 



We study the possibility that the light scalar mesons are q 2 q 2 states rather 
than qq. We perform a lattice QCD calculation of pseudoscalar meson scat- 
tering amplitudes, ignoring quark loops and quark annihilation, and find in- 
dications that for sufficiently heavy quarks there is a stable four-quark bound 
state with J PC = ++ and non-exotic flavor quantum numbers. 







I. INTRODUCTION 



The light scalar mesons have defied classification for decades |ffLH| . Some are narrow and 
have been firmly established since the 1960's. Others are so broad that their very existence is 
controversial. Scalar mesons are predicted to be chiral partners of the pseudoscalars like the 
pion, but their role in chiral dynamics remains obscure. Naive quark models interpret them 
as orbitally excited qq states. Others have suggested that they are q 2 q 2 @ or "molecular" 
states, [|J strongly coupled to tttc and KK thresholds. 

In this paper we propose a way to shed some light on the nature of the scalar mesons 
using lattice QCD. Previously scalar mesons have been treated like other mesons: their 
masses have been extracted from the large Euclidean time falloff of qq — qq correlation 
functions with the appropriate quantum numbers. Here we look for a ++ q 2 q 2 bound state. 
We construct q 2 q 2 sources, work in the quenched approximation, and discard qq annihilation 
diagrams so communication with qq and vacuum channels is forbidden. Also, we allow the 
quark masses to be large (hundreds of MeV), so the continuum threshold for the decay 
q 2 q 2 — > (qq)(qq) is artificially elevated. We then study the large Euclidean time falloff of a 
q 2 q 2 — q 2 q 2 correlator, looking for a falloff slower than 2rriq q , signalling a bound state. Such an 
object would have been missed by studies of qq correlators in the quenched approximation. 
We use shortcomings of lattice QCD to our advantage. By excluding processes that mix qq 
and q 2 q 2 , we can unambiguously assign a quark content to a state. The heavy quark mass 
suppresses relativistic effects, which we believe complicate the interpretation of light quark 
states. 

Our initial results are encouraging: within the limits of our computation we see signs of 
a bound state in the "non-exotic" q 2 q 2 channel, namely, the one with quantum numbers that 
could also characterize a qq state (1 = for 2 flavors, the 1 and 8 for 3 flavors). In contrast, 
the "exotic" flavor q 2 q 2 channel (1 = 2 for 2 flavors, the 27 for 3 flavors) shows no bound 
state. Instead it shows a negative scattering length, characteristic of a repulsive interaction. 
To obtain a definitive result will require larger lattices and more computer time, but this is 
well within the scope of existing facilities. 

In Sec. [TT]we give an overview of the ++ mesons. First we summarize the phenomenology. 
Then we summarize previous lattice calculations. We also review earlier studies of q 2 q 2 
sources on the lattice. [|5]|| Because these earlier works looked only at one (relatively small) 
lattice size they were unable to examine the possibility of a bound state. In Sec. |T| we 
summarize our computation. First we briefly review q 2 q 2 operators and discuss lattice size 
and quark mass dependence. Next, we review the improved lattice action we use to enable 
us to study larger lattices. Finally in Sec. |TV| we present our results and discuss their 
implications. We explore some of the directions in which our computation could be improved. 

A reader who wishes to skip the details can look immediately at Fig. |5| where we plot 
the dependence on lattice size of the binding energy associated with the exotic and non- 
exotic q 2 q 2 channels. The exotic channel shows a negative binding energy with the 1/L 3 
dependence expected from analysis of the (qq)(qq) continuum. H The coefficient of 1/L 3 
agrees roughly with Refs. [|5]|| and with the predictions of chiral perturbation theory. The 
non-exotic channel shows positive binding energy, but seems to depart from 1/L 3 , perhaps 
approaching a constant as L — > oo, which would indicate the existence of a bound q 2 q 2 state. 
Confirmation of this result will require further calculations on larger lattices. 
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II. OVERVIEW OF THE LIGHT SCALAR MESONS 



In this section we establish the context for our work. First we give a very brief intro- 
duction to the phenomenology of the lightest ++ mesons composed of light (u, d, and s) 
quarks. We give a sketch of the qq and q 2 q 2 models for ++ states and contrast them. More 
information can be found in Refs. [|l],[| and references quoted therein. Next we summarize 
existing lattice calculations which relate to the ++ channel. These fall into two classes: 
traditional searches for qq eigenstates and attempts to learn about low energy tttt scattering 
by studying q 2 q 2 sources. 

A. Phenomenology 

The known ++ mesons divide into effects near and below 1 GeV, which are unusual, 
and effects in the 1.3-1.5 GeV region which may be more conventional. Here we focus on 
the states below 1 GeV. Altogether, the objects below 1 GeV form an SU(3)f nonet: two 
isosinglets, an isotriplet and two strange isodoublets. The isotriplet and one isosinglet are 
narrow and well confirmed. The isodoublets and the other isosinglet are very broad and still 
controversial. 

The well established ++ mesons are the isosinglet /o(980) and the isotriplet ao(980). 
Both are relatively narrow: r[/ ] ~ 40 MeV, T[a ] ~ 50 MeVfj despite the presence of open 
channels (tttt for the fo and ttt] for the ao) for allowed s-wave decays. Both couple strongly 
to KK and lie so close to the KK threshold at 987 MeV that their shapes are strongly 
distorted by threshold effects. Interpretation of the fo and ao requires a coupled channel 
scattering analysis. The relevant channels are tttt and KK for the fo and ttt] and KK for 
the ao- In both cases the results favor an intrinsically broad state, strongly coupled to KK 
and weakly coupled to the other channel. The physical object appears narrow because the 
KK channel is closed over a significant portion of the object's width. No summary this brief 
does justice to the wealth of work and opinion in this complex situation. 

The other light scalar mesons are known as broad enhancements in very low energy 
s-wave meson-meson scattering. The enhancements are universally accepted, but their in- 
terpretation is more controversial. At the lowest energies only the tttt channel is open. The 
tttt s-wave can couple either to isospin zero or two. The 1 = 2 (e.g. tt + tt + ) channel shows 
a weak repulsion in rough agreement with the predictions of chiral low energy theorems. || 
The 1 = channel shows a strong attraction: the phase shift rises steadily from threshold 
to approximately tt /2 by ~ 800 MeV before effects associated with the fo complicate the 
picture. This low mass enhancement in the tttt s-wave is the a meson of nuclear physics and 
chiral dynamics. Recent studies support the existence of an S-matrix pole associated with 
this state at a mass around 600 MeV, which we will refer to as the cx(600). P,[10| The ttK 
s-wave is very similar to tttt. The exotic / = 3/2 (e.g. tt + K + ) channel shows weak repulsion. 
The non-exotic 1 = 1/2 channel shows relatively strong attraction. Black et al. identify 



We use the observed peak width into tttt and ttt] respectively, rather than some more model 
dependent method for extracting a width. 
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the enhancement with an S-matrix pole at approximately 900 MeV, which is known as the 
K(900).The enhancement is not in doubt, but the interpretation is, if anything, more con- 
troversial than the nn case.Q Other couplings of these objects (the a can couple to KK and 
the k can couple to rjK) are unknown because the relevant thresholds lie above the states. 
The large widths of these states reflect their strong coupling to the open decay channels nn 
and ttK respectively. 

The conventional quark model assigns the ++ mesons to the first orbitally excited mul- 
tiplet of qq states. As in positronium, ++ quantum numbers are made by coupling L = 1 
to S — 1 to give total J = 0. The ++ states should be very similar to the l +± and 2 ++ 
qq states that lie in the same family. These are very well known and form conventional me- 
son nonets (in SU(3){). Since they have a unit of excitation (orbital angular momentum), 
they are expected to be quite a bit heavier than the pseudoscalar and vector mesons. Most 
models put the qq ++ mesons along with their 2 ++ and 1 ++ brethren around 1.2-1.5 GeV. 

An idealized qq meson nonet has a characteristic pattern of masses and decay couplings. 
The vector mesons are best known, but the pattern is equally apparent in the 2 ++ or 1 ++ 
nonets. The isotriplet and the isosinglet composed of non-strange quarks are lightest and 
are roughly degenerate (e.g. the p and the uj). The strange isodoublets are heavier because 
they contain a single strange quark (e.g. the K*). The final isosinglet is heaviest because 
it contains an ss pair (e.g. the <fi). Decay patterns show selection rules which follow from 
this quark content. In particular, the lone isosinglet does not couple to non-strange mesons 
(<f) -/>■ 3tt). The mass pattern, quark content and natural decay couplings of a qq nonet 
are summarized in Fig. |l]a. These patterns seem to bear little resemblance to the masses 
and couplings of the light ++ mesons, a fact which led earlier workers to explore other 
interpretations. 

Four quarks (q 2 q 2 ) can couple to ++ without a unit of orbital excitation. Furthermore, 
the color and spin dependent interactions, which arise from one gluon exchange, favor states 
in which quarks and antiquarks are separately antisymmetrized in flavor. For quarks in 
3-flavor QCD the antisymmetric state is the flavor 3. Thus the energetically favored con- 
figuration for q 2 q 2 in flavor is (qq) 3 (qq) 3 , a flavor nonet. The lightest multiplet has spin 
0. Explicit studies in the MIT Bag Model indicated that the color-spin interaction could 
drive the q 2 q 2 ++ nonet down to very low energies: 600 to 1000 MeV depending on the 
strangeness content. || 

The most striking feature of a q 2 q 2 nonet in comparison with a qq nonet is an inverted 
mass spectrum (see Fig. [I|b). The crucial ingredient is the presence of a hidden ss pair in 
several states. The flavor content of (qq) 3 is {[ud], [us], [ds]}, where the brackets denote 
antisymmetry. When combined with (qq) 3 , four of the resulting states contain a hidden ss 
pair: the isotriplet and one of the isosinglets have quark content {udss, -^(uu — dd)ss, duss} 
and -^=(uu+dd)ss, and therefore lie at the top of the multiplet. The other isosinglet, uddu is 
the only state without strange quarks and therefore lies alone at the bottom of the multiplet. 
The strange isodoublets (usdd, etc.) should lie in between. In summary, one expects a 
degenerate isosinglet and isotriplet at the top of the multiplet and strongly coupled to KK, 
an isosinglet at the bottom, strongly coupled to 7T7t, and a strange isodoublet coupled to Ktt 



2 For example, the k(900) is not mentioned in Ref. [[l | 



3 



Mass 



ss 



KK 



us K* 



Ktt 



TTTTTT 



ud 
i 



1 



TTTT 



duss 



uu I dd 
ss 



udss 



KK 

irr) 



usdd 



Ktt 



udud 



J I I L 



TTTT 



1(a) 



1(b) 



FIG. 1. The mass pattern, quark content and natural decay couplings of (a) a qq nonet and 
(b) a q 2 q 2 nonet. 



in between (Fig. |l|b). The resemblance to the observed structure of the light ++ states is 
considerable. 

These qualitative considerations motivate a careful look at the classification of the scalar 
mesons. Models of QCD are not sophisticated enough to settle the question. For example, 
the q 2 q 2 picture does not distinguish between one extreme where the four quarks sit in the 
lowest orbital of some mean field, and the other, where the four quarks are correlated into 
two qq mesons which attract one another in the flavor {qq) 3 (qq) s channel. [11,3]] For years, 
phenomenologists have attempted to analyse meson-meson scattering data in ways which 
might distinguish between qq and q 2 q 2 assignments. A recent quantitative study favors the 
q 2 q 2 assignment. || However the qq assignment has strong advocates. [1^] We hope that a 
suitably constrained lattice calculation can aid in the eventual classification of these states. 



B. Existing Lattice Studies 



In this section we briefly summarize existing lattice calculations which bear on the clas- 
sification of the ++ mesons. There have been lattice studies of both the spectrum of ++ 
states and the mixing of qq states with glueballs. 

Unquenched spectroscopic calculations are just beginning to become available P~2| , |T3 . 
In principle, they are of interest because they would couple to a q 2 q 2 configuration if it is 
energetically favorable. One unquenched calculation reports tentative evidence of ++ state 
at an energy much lower than that reported in quenched calculations [|ll| . We return to this 
work briefly in our conclusions. Further insight from unquenched calculations will have to 
await more definitive studies. 

For the rest of this section, we restrict ourselves to consideration of quenched calculations. 
We will not discuss the mixing of glueballs with qq states, because we are interested in 
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FIG. 2. Quenched lattice calculations [15,17] of the ratio of the masses of various p-wave qq 
mesons to the vector meson mass, as a function of the pseudoscalar to vector mass ratio, ie with 
varying quark mass. The threshold for decay to P P is also shown. 



distinguishing qq from q 2 q 2 components of mesons. First we consider quenched studies of 
qq spectra. Then we describe attempts to extract meson-meson scattering lengths from 
quenched studies of q 2 q 2 sources. 



1. qq Spectrum Calculations 

The masses of the ++ qq states have been calculated on the lattice in the quenched 
approximation by various groups P^-4T7[. Some of their results are shown in Fig. As well 
as the J PC = ++ (a ), we have included data from the same groups on the other positive 
parity mesons J PC = l ++ (a 1 ) and 1 H Data for the 2 ++ (a 2 ) was not available. The 

spectra in Fig. [| behave roughly as qq states with orbital angular momentum should. In 
the heavy quark limit, as the pseudoscalar mass nip approaches the vector mass my, their 
masses approach one another because spin and spin-orbit splittings decrease with m q , and 
approach my because orbital excitation energy decreases with m q . To make this behavior 
manifest we plot rrtjpc /my versus mp/rriy, and note that nijpc /my approaches unity as 
mp/my increases. 



2. Pseudoscalar Scattering Length Calculations 



In the past, lattice studies of four-quark states have been undertaken in order to extract 
pseudoscalar-pseudoscalar (P-P) scattering lengths for comparison with the predictions of 



chiral dynamics. It is known || that the energy shift 5E of a two-particle state with quantum 
numbers a in a cubic box of size L is related to the threshold scattering amplitude, 



T 

-*- rv 
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SE. = E a - 2m P = ^ I 1 + 2.8373 ^ + 6.3752 )+■■■], (1) 

where mp is the mass of the scattering particles, and T a is the scattering amplitude at 
threshold in the channel labelled by a, which can be related to the scattering length, 

T a = -. 2 

m P 



For a more detailed discussion, see Ref. |18[]. In our case the channels of interest are exotic 
(I = 2, for two flavors) and non-exotic (I = 0, for two flavors). If the interaction is attractive 
enough to produce a bound state, then instead of eq. ([!]) one would find that SE goes to a 
negative constant as L — > oo. 

In order to distinguish between a bound state and the continuum behavior described by 
eq. (Jj), it is necessary to perform calculations for several different lattice sizes. Calculations 
with q 2 q 2 sources have been performed by Gupta et al. ||, who studied one lattice volume at 
one lattice spacing, and Fukugita et al. |§, who, for the heavy quark masses we are interested 
in, also studied only one lattice volume at one lattice spacing. Their results were therefore 
not sufficient to check the lattice-size dependence of energy of the two-pseudoscalar state, 
and investigate the possibility of a bound state. Our method follows theirs, but we have 
studied a range of lattice sizes. Their results are plotted along with ours in Fig. |. Where 
our calculations overlap, they agree. 



III. A q 2 q 2 EXERCISE ON THE LATTICE 
A. Quark contractions and flavor dependence 

For our purposes the salient categorization of q 2 q 2 correlators is into "exotic" channels 
(flavor states that are only possible for a q 2 q 2 state, 1 = 2 for two flavors, the 27 for three 
flavors) and non-exotic channels (flavor states that could be q 2 q 2 or qq, I = for two flavors, 
the 8 and 1 for three flavors). In the absence of quark annihilation diagrams, the 8 and 1 
are identical. When annihilation is included, the 1, like the 1 = for two flavors, can mix 
with pure glue. As shown in Fig. ||], the q 2 q 2 ++ correlation functions can be expressed 
in terms of a basis determined by the four ways of contracting the quark propagators ||18|| : 
direct (D), crossed (C), single annihilation (A), complete annihilation into glue (G). Since 
we are interested in q 2 q 2 states, we only study the D and C contributions. We will assume 
that all quarks are degenerate, so there is only one quark mass, and as far as color and 
spinor indices are concerned all quark propagators are the same. In our lattice calculation 
we will therefore build our q 2 q 2 correlators from color and spinor traces of contractions of 
four identical quark propagators, putting in the flavor properties by hand when we choose 
the relative weights of the different contractions. 

In the case of two flavors, there are two possible channels for a spatially symmetric 
source: 1 = 2 (exotic) and 1 = (non-exotic). Evaluation of the flavor dependence of the 
quark line contractions shows that the 1 = 2 channel is D — C, and I = is D + |C [18 . 
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FIG. 3. The four types of quark line contraction that contribute to the pseu- 

doscalar-pseudoscalar (P-P) correlation function. 

For three flavors, the possible channels are the symmetric parts of 3 x 3 x 3 x 3, namely 
1 + 8 (non-exotic) and 27 (exotic). As in the two-flavor case, the exotic channel is D — C. 
At sufficiently large Euclidean time separation, each contraction will behave as a sum of 
exponentials, corresponding to the states it overlaps with. Generically, all linear combina- 
tions will be dominated by the same state: the lightest. Only with correctly chosen relative 
weightings will the leading exponential cancel out, yielding a faster-dropping exponential 
corresponding to a more massive state. We will see in Sec. [TV] that the exotic (D — C) 
channel is the one where such a cancellation occurs, yielding a repulsive interaction between 
the pseudoscalars. For any other linear combination of D and C the correlator is therefore 
dominated by the lightest, attractive state. Without loss of generality, we can therefore 
study the following linear combinations: 

Exotic: J E = D-C 2 flavor: 1 = 2 3 flavor: 27 . . 

Non-exotic: J N = D + \C 2 flavor: 1 = 3 flavor: 1,8 

We conclude that if, as our results suggest, there is a bound q 2 q 2 state in the non-exotic 
channel, then this means that with two flavors, the 1 = channel is bound, and with 
three flavors both the 1 and 8 are bound. Once quark loops and annihilation diagrams are 
included, the 1 and 8 will split apart. Unquenched lattice calculations will be needed to see 
if they remain bound. 

B. Lattice action 

In our lattice calculations, we work in the quenched (valence) approximation, and use 
Symanzik-improved glue and quark actions. This means that irrelevant terms (0(a), 0(a 2 ), 
where a is the lattice spacing) have been added to the lattice action to compensate for 
discretization errors. 
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Improved actions are crucial to our ability to explore a range of physical volumes using 
limited computer resources. Because most of the finite-lattice-spacing errors have been 
removed, we can use coarse lattices, which have fewer sites and hence require much less 
computational effort: note that the number of floating-point-operations required even for a 
quenched lattice QCD calculation rises faster than a~ 4 . 

Improved actions have been studied extensively |7|JT9|-p2| , and it has been found that even 
on fairly coarse lattices (a up to 0.4 fm) good results can be obtained for hadron masses by 
estimating the coefficients of the improvement terms using tadpole-improved perturbation 
theory. For the energy differences that we measure, we find that the improved action works 
very well. There are no signs of lattice-spacing dependence at a up to 0.4 fm, so as well 
as greatly reducing the computer resources required, it enables us to dispense with the 
extrapolation in a that is usually needed to obtain continuum results. 

Our lattice glue and quark action parameters are summarized in table | and are described 
in detail in Ref. 0. For the glue we use a Liischer-Weisz (plaquette and 2x1 rectangle) 
action fl20| , pl~f . We measured the lattice spacing by NRQCD calculations of the charmonium 
P — S splitting, using the experimental spin- averaged value of 458 MeV. 

For the quarks we use a D234 action, which includes third and fourth derivative terms 
as well as an improved clover term. All the coefficients in the action are evaluated at tree 



level in tadpole-improved perturbation theory [23], using the mean link in Landau gauge to 



estimate the tadpole contribution. We work at a quark mass close to the physical strange 
quark: the pseudoscalar to vector meson mass ratio rapjray is 0.76. 

We collected data at two lattice spacings, a = 0.40 and 0.25 fm. The scaling of the 
hadron masses is good but not perfect: the pseudoscalar weighed 790 MeV on our coarser 
lattice and 840 MeV on the finer one. For our fits to eq. (ffj) in Sect. [IVl we used the average. 



C. Sources and fitting 

To look for bound ++ q 2 q 2 states, we investigate states of two pseudoscalar mesons on 
lattices of various volumes, keeping only quark-line-connected diagrams. We calculate the 
binding energy 8E^ in the exotic channel (flavor states that are only possible for a q 2 q 2 
state, ie I = 2 for two flavors, the 27 for three flavors) and the binding energy SE-^ in the 
non-exotic channel (flavor states that could be q 2 q 2 or qq, ie I = for two flavors, the 8 and 
1 for three flavors). We could have used sources based on preconceptions about maximally 
attractive channels in QCD. For example, one-gluon exchange and instanton interactions are 
known to favor the color 3 diquark channel, leading to interesting phenomenology at high 
density PH , but we verified that these have good overlap with the two pseudoscalar meson 
source that we used. 

Since we are interested in q 2 q 2 states, we only study the D and C contributions. For each 
gauge field configuration we evaluate the pseudoscalar correlator P(t) and the the direct 
("D") and crossed ("C") contributions to the two-pseudoscalar correlator. We use a wall 
source at t — 0, and both point and smeared sinks. 
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FIG. 4. Ratio of wall-source correlators R^s(t) (see eq. @), on a log scale, for lattice spacing 
a = 0.4 fm, lattice size L = 2 fm. The lattice was 40 spacings (10 fm) long, but only the first half 
is shown. For the smeared data we fitted to an exponential using t = 1 to 17. For the unsmeared 
we used t = 5 to 17. The fitted values of 5En agree. 



P(t) = Tr (G(t,x)G\t,i 

X 

D(t) = Y^[^(G(t,S)G^t, 

X 

C(t) = Tr (G{t,x)G*{t,x) G(t,x)G j (t, 



(4) 



where the trace is over color and spinor indices, and G(t, x) is the quark propagator in 
a given gauge background from a wall source at t = to the point x at time t. For 
smeared correlators, we performed covariant smearing at the sink. Note that the source for 
a pseudoscalar meson is ^75^, and the inverse propagator G _1 (t,x) = r y 5 G(t,x) , y 5 , so no 
factors of 75 appear in eq. (f|). 

We construct the exotic and non-exotic correlators 



Mt) 
Mt) 



(D(t)) 

W)) 



\{c{t)) 

(C(t)), 



(5) 



where the angle brackets signify an average over the ensemble of gauge field configurations. 

To obtain the binding energy 5E^ in the 1 = channel, and the binding energy SE-^ in 
the 1 = 2 channel, we construct ratios of correlators and fit them to an exponential 
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group 


P 


cr 1 (MeV) 


K 


mp/my 


m P (MeV) 


Gupta et al. || 


Wilson 6.0 


2590(60) |2| 


0.154 




940(30) 


Fukugita et al. || 


Wilson 5.7 


1220(180) m 


0.164 


0.740(8) 


620(90) 


this work 


L&W 1.719 


790(10) |] 




0.756(5) 


840(11) 




L&W 1.157 


495(4) 




0.756(4) 


790(6) 



TABLE I. Lattice parameters for studies of P-P scattering states. Lattice spacings are deter- 
mined by charmonium or upsilon P — S splitting 



(6) 



The ratios of correlators are expected to take the single exponential form only at large t, 
after contributions from excited states have died away. We followed the usual procedure of 
looking for a plateau, and checking that smeared and unsmeared correlators give consistent 
results. The results for a typical case are shown in Fig. f|. There is no difficulty in identifying 
the plateau and extracting 5E-^. 



IV. RESULTS AND DISCUSSION 
A. Our Results 

We measured 5En and SEe for several different lattice spacings and sizes. Our results 
are shown in Fig. [5] along with previous results from Refs. [|5],|6| . The exotic and non-exotic 
channels appear to scale differently as a function of L. The exotic channel falls like l/L 3 , 
which is the expected form for a scattering state, eq. ([I]). A fit is given in Table [TT], and shown 
in the figure. The non-exotic channel appears to depart from 1/L 3 falloff. To be complete, 
however, we have fitted the non-exotic data also to the form expected for a scattering state. 
The results are given in Table |I| and in the figure. 

The parameters of the lattice calculations are given in Table |. The lattice spacings were 
determined by quarkonium P—S splittings, using either charmonium or upsilon experimental 
measurements to set the scale. 

Although we only studied one quark mass, Gupta et al. H] repeated their calculation 
for a lower quark mass, corresponding to mp = 770 MeV, and found that SE-^ and 5E-^ 
were unchanged to within statistical errors. This suggests that studies of q 2 q 2 operators 
near qq-qq thresholds are not too sensitive to quark masses and leads us to combine the 
energy splittings from the different calculations in Table | on the same plot. For the 5Ee 
data we display both Gupta et al. and Fukugita et al.'s results with our own. We enlarged 
the error bar on Fukugita et al.'s point to 14%, since Ref. [^J quoted a 14% uncertainly 



in measuring their lattice spacing, which arises from the discretization errors involved in 
using an unimproved action on a coarse lattice. These are also apparent from the fact that 
Ref Pi's mp/my is close to ours, but its mp is significantly lower. For 5E^ we do not use 
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FIG. 5. P-P binding energy in non-exotic (N) and exotic (E) channels. The data at 
a = 0.08 fm are from ||; the data at a = 0.16 fm are from ||. The a = 0.25 fm and a = 0.4 fm 
points at L = 2 fm have been displaced slightly to either side in order to distinguish them. The 
lines are fits to eq. (see table ||). 
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channel 


amplitude T(MeV~ 2 ) 


X 2 /dof 


dof 


4/|T 


non- exotic 


1.20(5) x 1(T 5 


17.0 


5 




exotic 


1.18(3) x 10~ 5 


3.0 


6 


1.03(3) 



TABLE II. Threshold scattering amplitudes T obtained by fitting lattice calculations of 5 En 
and SEe to the scattering-state form eq. (|l|). In the chiral limit, 4/^T = 1. The non-exotic channel 
data does not fit the scattering-state form, so its fitted T is meaningless. 



Fukugita et al.'s data, since they included the annihilation diagrams, which we specifically 
exclude in order to see a q 2 q 2 state. 

Our results are consistent with those of Refs. |||§, even though we use much coarser 
lattices. This supports our use of Symanzik-improved glue and quark actions with tadpole- 
improved coefficients. As a further check on the validity of the improved actions, we note 
that at L = 2 fm, where we performed a calculation at two different lattice spacings for the 
same lattice volume, the results for the two lattice spacings agree very well. There is no 
evidence of any discretization errors. 

For the exotic q 2 q 2 system, the fit to eq. (||) is quite good, and the fitted scattering 
amplitude is remarkably similar to the result expected in the chiral limit, 4/J.T = 1.0 We 
conclude that there are no surprises in the exotic channel - the interaction near threshold 
appears repulsive and the strength is close to that predicted by chiral perturbation theory. 

The non-exotic q 2 q 2 system, however, does not fit the expected scaling law at large 
L. The fit to eq. ([]]) has a very large x 2 ? an d is so poor that the extracted amplitude 
T is meaningless. Instead 5E^ appears to be approaching a negative constant at large L. 
Instead of a scattering state, we appear to be seeing a bound state in the non-exotic channel. 
Although our data are suggestive, they are not conclusive. It would be very interesting to 
gather more data at L > 4 fm, as well as at a range of quark masses, in order to verify the 
existence of this new state in the quenched hadron spectrum. 

B. Interpretation and Discussion 

We have found evidence for a q 2 q 2 bound state just below threshold in the non-exotic 
pseudoscalar-pseudoscalar s-wave. In 2-flavor QCD the bound state would correspond to 
an isosinglet meson coupling to ttti. In 3-flavor QCD the non-exotic channel corresponds 
to an entire nonet including two non-strange isosinglets and an isotriplet, and two strange 
isodoublets (see Fig. [I|b). We work with a large quark mass so our results are not directly 
applicable to tttt scattering, but they do resemble physical KK scattering.^] The known 
isosinglet /o(980) and isotriplet a (980) mesons are obvious candidates to identify with the 
non-exotic q 2 q 2 bound states we seem to have found on the lattice. 



3 Since we did not calculate fp at our quark masses, we have used the value fp = 148 MeV, 
derived from Ref. M, Table 1. 

4 Although we work in the SU(3)f limit where all quark masses are equal. 
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FIG. 6. Binding energy in MeV of the two-particle state in our toy model eq. (|^), as a function 
of particle mass in GeV. 

We believe the quark mass dependence of the non-exotic q 2 q 2 state is quite different 
from a standard qq lattice state. In the quenched approximation the masses of qq states like 
those shown in Fig. [2] are roughly independent of mp. Note especially that the masses are 
smooth as they cross the threshold, 2rrip. In contrast, we believe that the q 2 q 2 state we may 
have identified is strongly correlated with the PP threshold when the quark mass is large, 
and departs from it in a characteristic way as the quark mass is reduced. (Indirect support 
for this comes from Gupta et al.'s finding that their binding energy is independent of the 
pseudoscalar mass.) In particular, we believe that the bound state will move off into the 
meson-meson continuum as mp is reduced toward the physical pion mass. 

To explore the mp dependence of our results, we have made a toy model based on a 
relativistic generalization of potential scattering. We write a Klein-Gordon equation for the 
s-wave relative meson- meson wavefunction, 0(r), 

- 0"(r) + (2m P - U{r)) 2 4>{r) = £ 2 0(r), (7) 

with the boundary condition that 0(0) = 0. For U(r) = the spectrum is a continuum 
beginning at E — Imp as required. In the non-relativistic limit mp <C |Z7|, eq. (0) reduces 
to the Schrodinger equation with an attractive potential —U(r) (for U(r) > 0). For sufficient 
depth and range, this potential will have a bound state. However, as mp — > 0, the potential 
term in eq. (0) turns repulsive and the bound state disappears. Thus, if one keeps the depth 
and range of U fixed as one decreases mp, the bound state moves out into the continuum 
and disappears. To be quantitative, we have taken a square well, U(r) = Uq, for r < b, and 
U(r) = for r > b. We chose a range b = l/m n w 1.4 fm, and adjusted Uq such that the 
bound state has binding energy of 10 MeV when mp ~ 800 MeV. The bound state does 
indeed move off into the continuum (first as a virtual state) when mp goes below 330 MeV. 
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The behavior of the bound state in this toy model is shown in Fig. || Note this toy model 
is not meant to be definitive^ but it illustrates the expected behavior of a P P bound state: 
tracking 2mp with roughly constant binding energy as mp falls, then unbinding at some 
critical mp. 

On the basis of our lattice computation and the mp dependence suggested by our toy 
model, we believe it is possible that all the phenomena associated with the light scalar 
mesons are linked to q 2 q 2 states. The narrow ++ isosinglet / (980) and isotriplet ao(980) 
mesons near KK threshold can be directly identified with q 2 q 2 lattice bound states (top line 
of Fig. [TJb). The broad k(900) and <r(600) (middle and bottom lines of Fig. |I|b) couple to low 
mass (nn or irK) channels. We speculate that they are to be identified as the continuum 
relics of the same objects which appear as bound states of heavy quarks. 

Of course, a thorough examination of this question would require implementing flavor 
SU(3) violation by giving the strange quark a larger mass. This would mix and split the 
isoscalars, shift the other multiplets (see Fig. |I|b), and dramatically alter thresholds. For 
example, the 1=1 q 2 q 2 state couples both to KK and nr] (through the ss component of 
the 77) in the quenched approximation. The fact that the physical KK and irrj thresholds 
are significantly different would certainly affect the manifestation of bound states such as 
those we have been discussing in the S[/(3)-flavor-symmetric limit. 



C. Conclusions and Future Work 

We have presented evidence for previously unknown pseudoscalar meson bound states in 
lattice QCD. Our results need confirmation. Calculations on larger lattices are needed, and 
variation with quark mass, lattice spacing, and discretization scheme should be explored. 

In the real world a ++ q 2 q 2 state may, depending on its flavor quantum numbers, mix 
with ++ qq and glueball states. It seems natural to expect that for sufficiently heavy quarks 
a bound state will remain, but only full, unquenched lattice calculations can confirm this. 

It is possible that an existing unquenched study of ++ qq operators may show some 
corroboration of our results. In Ref. [|b|] the authors study qq sources with dynamical 
fermions. Although their interest was in exploring the mixing of ++ qq with glueballs, 
there is nothing to stop their qq source from mixing with q 2 q 2 . So they should be sensitive 
to the q 2 q 2 bound state we have identified. It is therefore quite interesting that they report 
a ++ state with an anomalously low mass ~ 800 MeV. 

If light q 2 q 2 states are, in fact, a universal phenomenon, and if the cx(600) is predominant- 
ly a q 2 q 2 object, then the chiral transformation properties of the a have to be re-examined. 
The 7r and the cr(600) are usually viewed as members of a (broken) chiral multiplet. In the 
naive qq model both n and a are in the (~, |) © (|, \) representation of SU(2)l ® SU{2)p 
before symmetry breaking. In a q 2 q 2 model, as in the real world, the chiral transformation 
properties of the a are not clear. 



5 We could have chosen a different relativistic generalization of the Schrodinger equation which 
would have preserved the bound state as mp — > 0. For example, we could have replaced (2mp — U) 2 
by 2m p — 2mpU\ — ?7f , and fine-tuned U\ and U2 to provide binding at arbitrarily low mp. 
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If the phenomena that we have discussed survive the introduction of differing quark 
masses, then they will also have implications for heavy quark physics. For example, there 
could be a ++ bound state just below the decay threshold for two D mesons in the char- 
monium spectrum. 

Finally, we note that calculations similar to ours could be undertaken in the meson- 
baryon sector and in other J PC meson channels. It has long been speculated that the 
A(1405) is some sort of KN bound state and q 2 q 2 states have been postulated in other 
meson-meson partial waves. 
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